Ratcheting up energy by means of measurement 
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The destruction of quantum coherence can pump energy into a system. For our 
examples this is paradoxical since the destroyed correlations are ordinarily considered 
!>■ , negligible. Mathematically the explanation is straightforward and physically one can 

identify the degrees of freedom supplying this energy. Nevertheless, the energy input 
can be calculated without specific reference to those degrees of freedom. 
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Introduction. 



Under some circumstances, a measurement can put energy into a system. This has long been 



\ recognized, for example, in locating a quantum particle "under" a barrier. However, for the 
process described below, the decohering of separated systems, one might have thought it to 
be innocuous, something that takes place without external intervention. We show that in 
fact the energy increases due to this decohering. Moreover, this has a classical analogue, so 
^ I that one can say that the energy influx arises because of compliment arity (in a sense to be 
■ explained below). 

The system is a pair of oscillators, A and B, and could represent an atom and a field or 
two sorts of oscillators whose self-interactions are substantially harmonic and which couple 
linearly when close. The excitations are called bosons. We take as Hamiltonian 

H = Ho + VsB = ujatt^a + cobb^b + g {a^ + a) {b^ + b) , (1) 

in the usual notation. This is the spin boson model [1]. As we proceed we will also refer to 
recent results on the Jaynes-Cummings model [2], i.e., with coupling Vjc = g {a^b + b'^a). 

We study the coming together of the two systems for a brief time, followed by their 
separation (for example, a pulse of light impinging on an atom). Thus one could consider 
(7 to be a function of time. Initially we take the density matrix to be a product, p(0) = 
Pa(0) ® Pb(0); the systems are not entangled. Subsequent to the encounter the state is in 
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general entangled; however, once they are separated the correlations associated with the 
entanglement can be dropped (provided one does not do an EPR experiment) and the true 
time-evolved p(t) = exp{—iHt)p{0) exp{iHt) can be replaced by effective, individual density 
matrices, Pa(t) = Tib pit) and pb(t) = TTap{t). 

For both Vsb and Vjc, diagonalization of H is straightforward. There is however a 
significant difference that is reflected in the temporal evolution of the operators as well as 
in the density matrices. For Vjc there is boson conservation, which implies: 1) The images 
of and under time evolution are linear combinations of the time-0 quantities (a and b 
do not enter). 2) If Pa{0) and ^^(0) are diagonal in the number operator basis, they stay 
that way. 3) On the other hand, coherent states, which are of the form exp(^aa''' + Zbb'^)\0), 
are mapped into other coherent states (there is a unitary transformation on the z's) and 
thus preserve non-entanglement. Now for Vsb there is also a rotation that diagonalizes H; 
after all, as oscillators the coupling is of the form ^XaXb (where a = [xasf^^iVaj \f^a)l ^f^-, 
etc., and 7 = const). However, because of the different a;'s that enter for A and i3, the 
time-transformed is a linear combination of all four operators, a\ a, b^ and b. It follows 
that the non-entanglement of coherent states is lost, Pa{t) and pb{t) are not diagonal in the 
number (Fock) representation, and the number of bosons is not conserved. 

We next consider a more elaborate situation. Imagine a collection of ^'s and H's that 
repeatedly come in contact. Thus Ai and Bi meet and separate, subsequent to which Ai 
goes on to encounter a different B, say B2. The usual way to treat the second encounter is 
to use, for the state of ^1, the reduced density matrix, Pa{t), from its last encounter. If B2 
also emerged from a similar encounter our estimate of its state is Pbit)- So in principle we 
should maintain a collection of Pa's and p^'s, and, perhaps randomly, allow pairs to interact. 
In practice, this procedure does not affect our conclusions and we use the simpler method 
of taking the output of one encounter, both Pa{t) and Pb{t), and using it as the input for the 
next. Thus the initial density matrix for the next encounter is not p{t), but Pa{t) ® Pbit)- 
Based on the assumed multitude of .4's and B's and the unlikehhood of a pair immediately 
re-encountering one another, this replacement should make no difference. 

The surprising result is that it makes a big difference. We will explain this both mathe- 
matically and physically, but we first present the results. Fig. 1 shows the probability dis- 
tributions of boson number after 25 successive encounters, with parameters cUa — i, — 
g = 0.2 [3], t = 4, initial states, — 2, rib — 1) and cutoff 20 (i.e., 21 oscillator levels are 
allowed). The long-time distributions are exponentials with A and B having roughly the 
same dropoff behavior, implying that distinctly different amounts of energy are sequestered 
in the two modes. For some runs, however, there is a slow drift in the A and B values, so 
that this may be a transient (and numerical considerations prevented our checking). The 
equality of (ua) and (n^) was observed also for the Jaynes Cummings model [2], and to much 
greater precision (without the drift that makes us here suspect transience). The remarkable 
feature of Vsb, not observed for Vjc, is that the average excitation of both oscillators grows 
as a function of encounter, which is to say, as a function of time. In Fig. 2a is shown the 
bosonic content for the same parameters as in Fig. 1. In Fig. 2b is a run with different 
parameters [ujb = 3, g = 0.5, t = 15). In this case A and B seem not to tend to the same 
limit although both increase. This growth is apparently linear and certainly represents an 
increase in the total energy in the system. Although one usually thinks of the destruction 
of quantum correlations as purely an information issue, here it has direct consequences. 

Mathematical explanation. 

Mathematically, we replace a distribution of two variables by its marginals; it's more compli- 
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FIG. 1: Probability distribution of bosons after 25 "encounters" for systems A and B (marked pa 
and pb). Note that they are essentiahy the same. 
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FIG. 2: Value of the average excitation level in system A (o's) and B (x's). Fig. 2a uses the 
parameters of Fig. 1. For Fig. 2b shghtly different parameters are used, with longer contact times 
and stronger coupling. In this case (ua) and (n^) grow, but do not tend to a common value. 

cated than in classical probability, since the off-diagonal elements of p are complex, but the 
principle is the same. How this causes trouble can be seen almost without calculation. Let 
^(0) = \na,nb) = A/'at""6t"''|o), with M = I/Vnjn^. Under time evolution the operators 
become linear combinations of a^, a, etc. Thus 

= J\f [ttia^ + Q!2a + a^b^ + 046]'*" 

X [Pia^ + 132a + /?3&^ + /?4&] |0) , (2) 

where ak — o-kit) and = A; = 1, . . . , 4, are explicit functions of time [4]. Remarkably, 

Eq. (2) implies that under the exact dynamics this state never has more than Ua + Ub bosons 
of either type. Thinking perturbatively it would seem that more bosons could be created, 
but this is not so. However, if one takes ilj{t)ilj^{t) (using Eq. (2)) and forms from it pa and 
Pb, then the correlations of the exact dynamics are lost, and the perturbative considerations 
apply. Thus the term b^ + b in the Hamiltonian is applied to pb irrespective of what happened 
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to A. This allows unlimited numbers of bosons to be created. (This also lowers n, but there 
is an asymmetry in that n > 0.) 

A second way to see the increase is through short-time behavior. We give an abstract 
statement of the problem, incidentally showing the phenomenon to be more general. Let 

H = Ha + Hb + Va^Vb, (3) 

where A and B need not be oscillators. At t — the density matrix is taken to be p(0) = 
Pa{0) ® Pb(0). Under time evolution p{t) = exp{—iHt)p{0) exp{iHt), which is in general 
entangled; pa(^) = Tr^ p{t) and similarly for B. Define Ap{t) = pA{t) ® — p{t). It is 
easy to show that 

AH = Tr [Ap{t)H] = IV [Ap{t) {Va ® Vb)] , (4) 

which is to say that the expectation of the free Hamiltonian, Ha + Hb, is not affected by 
the replacement. Thus it is the interaction term that gives rise to the effect exhibited in Fig. 
2. For the next step we calculate the short time form of p(t), compute its marginals and 
evaluate the difference. After some calculation (which we will present in [5]), one obtains 

A// = ^{ [Tr^ [VIpa{^)) - (Tr^ (VaPa(O)))'] 

^TiB{[VBAHB.VB\]pBm 
+ [Tib (Vjpi?(0)) - (Tr^ (V^bPb(O)))'] 

xTYA{[VAAHA,VA]]pAm] 
+0{e) . (5) 

Note that [Tr^ {VIpa{^)) - (Tr^ (KiPa(0)))^] > 0, and can only be zero if pA is concentrated 
on a single value of Va- Therefore to show AH > we examine the commutators. 

The example of interest is when A and B are harmonic oscillators. We use a coordinate 
representation. Thus Ha = Pa/2 + ^\x\/2, Va = \/\g\xA, etc. By direct calculation one 
finds [Va, [Ha-, Va]] = Ififl > [6]. This shows that the replacement of p by pa®P_b necessarily 
increases the expected value of the Hamiltonian, for short times. 

Physical explanation. 

The paradoxical aspect of our result is that one expects that in a gas of interacting particles 
there is little physical significance to their, say, momentum correlations, once they have 
separated a significant distance. This should apply even more to quantum coherence. It is 
true that there is a continuing loss of information, but losing quantum correlations should 
not heat the gas. You don't burn your finger because of a partial trace over a density matrix. 

Two phenomena shed light on this situation. First there are the results of entropy pro- 
duction in computing, namely the fact that the one can avoid any thermodynamic cost in 
a computation, provided one does not erase [7]. We will see how this plays a role. Also 
relevant are physical models of ratchets [8-14], the relation to which will disabuse anyone of 
the idea that this system could give rise to a perpetuum mobile. 

As remarked, the coupling, g, can be thought of as a function of time. In fact, it must 
be, since if the coordinates Xa and a;^ arc physical coordinates then true oscillators will 
continue to interact at all distances — the farther, the stronger. Alternatively, one could 
think of these oscillators as internal coordinates on particles with physical position r{t)- 
Then the physical approach and separation of the particles leads to a couphng coefficient of 
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the form g{ra{t),rb{t)), with g ^ as \ ra{t) — rb{t)\ grows. Similar considerations apply 
if the "oscillator" is a mode of a field, although sometimes one can make this idealization 
without running into trouble [2]. 

Nevertheless, in the situations contemplated here the time-dependence of the coupling 
constant implies that energy conservation need not apply. For this reason we recall the 
concept of ratchet, where the turning on and off of a potential induces directional flow. In 
biological applications this requires external energy, which is the role of ATP. In our case, a 
full explication will depend on the physical system that Eq. (1) represents. Suppose that the 
oscillator is an internal particle coordinate, borne by the translational degrees of freedom 
of that particle from place to place, repeatedly encountering other particles. Then as two 
oscillator-bearing particles approach one another, the coupling energy g{ra{t),rb{t))xaXb 
begins to affect the translational motion of the particles themselves. Since the entire system 
can be described by a time-independent Hamiltonian, the source of the energy that enters 
the oscillator coordinates is necessarily the translational degrees of freedom of the moving 
particles. If too much is withdrawn, the particles will cease encountering one another (or 
even separating), and the tracing over the "other" degree of freedom inappropriate. 

The fact that the passage of energy to internal degrees of freedom can cool a gas (trans- 
lationally) is no surprise; what is of interest is that it comes about through the destruction 
of information, in this case the destruction of quantum coherence. This kind of information 
loss is not ordinarily considered a source of energy transfer, although, as alluded to above, 
erasure can have thermal consequences. 

Another argument highlighting the paradoxical nature of our result is that we calculate 
the energy increase without any explicit model of particle interactions. For the realization 
mentioned above, g{ra{t),rh{t))XaXb, one would have expected that the dependence of g on 
distance should play a role in understanding the energy transfer, but somehow we don't 
need that; decoherence alone drives the process. 

Classical oscillators. 

To broaden the perspective we mention a classical analogue. Imagine two oscillators with 
frequencies cua and cub- At a random time couple them by adding jXaXb to the Hamiltonian 
and allow them to evolve, again for a random time. Repeat the coupling and uncoupling 
many times. The result is that the energy grows exponentially. This phenomenon is easy to 
check numerically, but it is also easy to develop a qualitative explanation. The system moves 
on a 2-dimensional torus in 4-dimensional phase space; at any moment it is on a particular 
torus. When the coupling is switched (adding or subtracting ^yXaXb), the collection of tori 
changes and the system point continues on the torus associated with the new dynamics. The 
new torus need not have the same energy as the old one. The motion can be described as 
a jumping from torus to torus, and as a consequence, a random walk in energy. However, 
because the system is strictly linear, the dynamics is scale-invariant. Therefore the natural 
description of the random walk is not in energy, but in its logarithm. That a drift-free random 
walk in logE leads to exponential growth in E can be seen as follows. Let u = logE and 
suppose that the distribution function for u is p{u) = exp ((-u — /2a'^) /a\/2Tx. For this 
p, (u) = Uq. But the expectation of E is not exp(Mo); rather {E}/ cxp {{logE)) = exp(cr^/2). 
Under diffusion for time t, with diffusion coefficient D, the spread, a^, of the distribution 
changes to o"^ + 2Dt. Since the expectation of the central point Uq does not change, the 
expectation of E grows like exp{Dt). As a simpler example on can imagine a single oscillator 
that switches at random times between two frequencies cu and uj'. Now the "jumping" is 
between two ellipses in phase space and the energy grows exponentially. 
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This classical analogue suggests that one can think of the quantum phenomenon as a 
manifestation of complimentarity. For the quantum system it is not possible to retain all 
information after a measurement, so that restarting the system with partial information 
causes it to resemble a classical system for which correlation information has been destroyed 
"by hand." (In the work just described we accomplish this by letting the system run for 
"random" times. With non-random times it will sometimes grow exponentially, sometimes 
oscillate.) 
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